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lent agreement with the numerical solutions especially for the case where the 
initial phase difference between the bubbles is small. With a step-function 
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Abelian field strength and other gauge-invariant quantities are obtained in 
closed form. Possible extensions of the solution to the case of the electroweak 
phase transition and the generation of primordial magnetic fields are briefly 
discussed. 



°Electronic address: e . j . copelandQsussex .ac.uk 
b Electronic address: p.m.saffin@damtp.cam.ac.uk 
c Electronic address: o.tornkvist@damtp.cam.ac.uk 



1 



1 Introduction 



Phase transitions are a universal feature in field theories. For example, they may explain 
how quark confinement has emerged and possibly how the electroweak symmetry was 
broken. It is also through such transitions that topological defects form in condensed- 
matter systems in the laboratory and, we believe, could have formed in the early universe 
(see 0] and [|| for reviews). Phase transitions of the first order play a particularly 
significant role in particle cosmology, either as a progenitor of inflation |3|, |J |5j or as a 
source of the baryon asymmetry associated with the electroweak transition ||. Popular 
models of electroweak physics such as the Minimal Supersymmetric Standard Model 
exhibit a first-order transition for a range of experimentally allowed parameters . First- 
order phase transitions are driven by the nucleation of bubbles of true vacuum within 
a sea of false vacuum. The bubbles, responding to the pressure caused by the energy 
difference of the two vacua, expand and subsequently collide. Over the years, a number 
of authors have turned their attention to the dynamics of such bubble collisions in gauge 
theories || ^ [10], [11], 0, J3|. In spite of this effort, it is fair to say that much remains 



unclear about the role the various fields play in the collision process. This should not 
be surprising, as we are here dealing with a non-linear system which is difficult to treat 
analytically. 

Much of the interest in bubble collisions derives from the pioneering work of Kibble 
H , who demonstrated the possible formation of topological defects in a cosmological 
phase transition for cases where the vacuum manifold of the Higgs potential is non-trivial. 
In particular, defects such as strings (or monopoles) may be produced in the collision 
of three (four) bubbles, provided that inside each bubble the scalar Higgs field takes a 
different random value on the vacuum manifold. In Ref. [|IJ[] the mechanism for defect 
formation was formulated in terms of a geodesic rule, according to which the phase (or, 
more generally, the isospace orientation) of the scalar field between two colliding bubbles 
takes values along the shortest path on the vacuum manifold that connects the values 
in each bubble. The geodesic rule has since been verified for defect formation in various 



models with broken global symmetry [15, 16|. In gauge theories, however, Rudaz and 



Srivastava || have argued that the geodesic rule does not apply, since the phase of the 
scalar field has no gauge-invariant meaning. Instead, they suggested that the dynamics of 
the fields would play an important role in the production of gauge defects in a first-order 
phase transition and could possibly lead to a suppression in the rate of defect formation. 
This claim was subsequently investigated by Hindmarsh, Davis, and Brandenberger 



T0|| . Assuming the gauge fields to be negligible at the moment of collision, they obtained 



solutions which appeared to adhere to the prescription of the (naive) geodesic rule and 
found no suppression in the probability of defect formation. Later, however, Copeland 
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and Saffin [Iz] demonstrated unequivocally that the geodesic rule can be violated in a 
U(l) gauge theory as a consequence of the detailed energetics in bubble collisions. The 
reason is that violent fluctuations in the modulus of the scalar field cause symmetries 
to be restored locally, allowing the phase to "slip" by an integer multiple of 2ir and 
thus change the number of string defects. In fact, when the initial phase difference of 
the bubbles and their initial separation are sufficiently large, string defects can arise 
when only two bubbles collide, unlike the usual picture of string formation obeying the 
geodesic rule, which would require at least three bubbles. In this article, however, we will 
be dealing with phase differences and bubble separations which are too small to cause 
defects to form through violation of the geodesic rule. 

An appealing aspect of first-order transitions involving gauge fields is that the bubble 
collisions give rise to magnetic fields. These are produced by a gauge-invariant current 
which develops across the surface of intersection of the two bubbles due to the gradient of 



the phase of the scalar field ||11|| . Because bubble collisions occur out of thermal equilib- 
rium and the bubble walls possess large kinetic energy, the magnetic fields thus created 
can be considerably larger than the thermal background of magnetic-field fluctuations. 
If the energetics permits the geodesic rule to be violated, these magnetic fields become 
associated with the formation of local (gauged) defects. 

Magnetic fields could also be generated in extensions of the SU(2)xU(l) Standard 
Model with a first-order phase transition and have been studied by a number of authors 



17], [T8|, |l9fl. Such fields, created in the early universe, could evolve into seed fields for a 



galactic dynamo mechanism that would in turn amplify them to the strength observed 
in galaxies today. The common view has been that it is generally too difficult to produce 
large enough fields at the electroweak phase transition on sufficiently large length scales, 
but recently Davis, Lilley and Tornkvist |2(| have argued that in the context of a universe 
with less matter than the critical density (and particularly with a positive cosmological 
constant) a much smaller primordial field would be sufficient for explaining the present 
galactic fields. The reason is simple: The increased age of galaxies in such a universe 
gives the primordial seed field more time to grow in strength. 

In this article we shall concentrate on the case of the U(l) Higgs model, partly to 
simplify the analysis but also to compare our results with those of Kibble and Vilenkin 
|TT|] , who developed a formalism for calculating the magnetic fields generated in bubble 



collisions. Of course, one should not really call it a magnetic field in the case of an 
Abelian gauge theory, owing to the lack of a massless photon in the broken-symmetry 



phase. Instead we will generally refer to it as a "field strength". The authors of Ref. [ 1 1 



made some potentially severe approximations in order to solve the field equations: The 
modulus of the complex Higgs field was kept fixed in the true vacuum, and the initial 
phase of the Higgs field at the moment of collision was approximated by a step function. 
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With these basic assumptions, they were able to derive some fascinating analytical results 
concerning the generation of the field strength in the collision region. They demonstrated 
how the gauge-invariant current across the interface of the two colliding bubbles produces 
a ring of flux of the Abelian field strength which encircles the bubble intersection region. 
On top of this they defined a gauge-invariant phase difference between the two bubbles 
and solved for its future evolution to show that it oscillates and eventually relaxes to zero 
as the bubbles merge. This example of phase equilibration between two colliding bubbles 



is in fact a generic feature of models with either local or global symmetry |Tl], |21 



Here we go beyond the approximations introduced in Ref. JTIJ and, by doing so, test 
their regime of validity. In particular, we allow for a smoother interpolation of the phase 
between the two bubbles, using an approximation to the bounce solution for a single 
bubble, and the fact that we can represent the scalar field of two well separated bubbles 
as a linear superposition of the scalar fields for each bubble [10, 13]. Remarkably, we can 
obtain explicit time-dependent expressions for both the gauge-invariant phase and the 
Abelian field strength. We compare these expressions to those in Ref. JI1J as well as to 
the numerical solution of the field equations. The results are very encouraging; we find 
incredible agreement between our analytical solutions and the numerical solutions. We 
also demonstrate to what extent fields in the collision region are well approximated by 



the step-function ansatz adopted in Ref. [11 



Real phase transitions naturally possess lots of complicating factors which we do not 
take into account here, such as frictional damping of the bubble-wall velocity. Particles 
that acquire a mass due to couplings to the Higgs field are massless outside the bubble 
and massive inside. Not all the particles external to the bubble have enough energy to 



become massive and must reflect off the bubble wall, slowing it down [£lj, ^2 



Another complication is the conductivity associated with the plasma surrounding the 
expanding bubble wall. Charged particles in the plasma couple to the gauge fields and 
cause them to lose their energy to the thermal bath. This has the effect of reducing the 
likelihood of forming loops of flux during the bubble collision as pointed out in [JlO, [11 
and confirmed in | 12| . Furthermore, the existence of a highly conductive plasma has the 
effect of conserving magnetic flux by virtue of Faraday's and Ohm's laws. 

In spite of these potential problems, we believe that the simple system of colliding 
bubbles is well worth investigating, not least because the field theory is highly non-linear 
and yet we are able to find extremely accurate analytical solutions. Moreover, as we 
shall discuss in the conclusions, the method and results can be generalised and applied 
to different areas, such as the Standard Model and its extensions. 

The paper is organised as follows. In Sec. 2 we review the standard results due to 
Coleman [23 on bounce solutions and bubble nucleation in first-order phase transitions. 
Following Ref. |25[ , we also investigate the evolution of bubbles of arbitrary nucleation 



4 



radius and find that all expanding bubbles are well approximated by the bounce solution 
after the bubble size has increased only by a small amount. In Sec. 3 we review the 
key approximations and results obtained in Ref. [IT], in particular the solutions for the 



gauge-invariant phase and the field strength arising from the step- function approximation. 
We rewrite these solutions, and solutions for other gauge-invariant quantities, as closed- 
form expressions in terms of Bessel functions. In Sec. 4 we obtain improved solutions 
by using the smooth bounce solution for the scalar field of each individual bubble and 
the linear superposition ansatz for the combined scalar field. We derive explicit time- 
dependent solutions for the gauge-invariant phase and the Abelian field strength. The 
numerical lattice computation is described in Sec. 5. In Sec. 6 we compare the numerical 
field evolution to our new smooth solution and to the results obtained in Ref. [11]. Of 
particular note is that the step-function approximation leads to discontinuities in the field 
strength and in the phase that are not there in the smooth approximation. Moreover, the 
energy density associated with the gauge field arising from the step-function approach 
is about sixteen times larger than the true solution found both numerically and in the 
smooth approximation. We conclude in Sec. 7. 

2 Kinematics of Expanding Spherical Bubbles 



In a pioneering paper Coleman [M] investigated the physics of bubble nucleation in the 
decay of the metastable vacuum. He obtained the so-called bounce solution for the mod- 
ulus of the scalar field corresponding to a single nucleated bubble. If the bubble radius is 
large compared to the wall thickness, the bubble wall can equivalently be described as a 
relativistic membrane. The motion of such a membrane subjected to the outward pres- 
sure caused by the energy difference of the two vacua can be derived and was investigated 



in the spherical case by Garriga |[25|| . Below we will reproduce a number of key results 
in Refs. [p4] , ^5] , extracting details of the bubble- wall dynamics that we shall require 
later. Of particular use will be an analytic approximation, Eq. (|]), to the actual bounce 
solution, Eq. (|), as well as the most probable size of the nucleated bubble, Eq. (|10D, and 
the associated wall thickness, Eq. (|i~l|). 

We shall see below that bubbles created via an instanton description of tunnelling have 
0(1,3) symmetry. When there are two bubbles present the spatial rotation symmetry is 
reduced to 0(2), so in fact the two-bubble collision has an 0(1,2) symmetry ||. When 
we have two bubbles there is the chance that they may nucleate at different times and so 
have different sizes when they collide. We may however always boost to a frame where 
they nucleate simultaneously; this is the frame in which we shall work. 
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2.1 Bubble Nucleation 



For simplicity we choose a Higgs potential^ that is cubic in 



V(\<>\) = M\^-r ] ri^ + e)-e\r ] \ (1) 




where is the scalar Higgs field and < e < |. It has a metastable minimum at = 
and a degenerate global minimum at = rje 19 , 9 e [0,27r), separated by a barrier of 
height [4(1 + e) 3 /27 — e] A?] 4 . The difference in energy density between the two minima 
is p = eXr] 4 , and this constant has been subtracted in Eq. ([]]) in order to obtain zero 
energy density in the metastable minimum. The Higgs-boson mass Mh in this model is 
given by M\ = 4(1 + e)Xr) 2 . 

The nucleation of a bubble of true Higgs vacuum \<f>\ — r\ occurs by means of tunnelling 



from = through the barrier [21]. The tunnelling is expected to be dominated by so- 
lutions to the Euclidean equations of motion, with a nucleation probability oc exp(— S E ), 
where Se is the Euclidean action. S E is given by 



S E [0,AJ = / d A x 



(2) 



where a positive definite metric is used. The covariant derivative is = <9 M + ieA^ 
with Ffj, v = d fl A u — duA^, and A^ is the associated U(l) vector potential. Tunnelling is 
therefore dominated by those instantons with lowest Euclidean action, something which 
is achieved by having vanishing vector potential and a constant phase for the scalar field 
0. 

In the absence of gauge fields the instantons which have lowest action (0) are those 



with 0(4) symmetry, i.e. = <p(s), where s — yt| + x 2 and i E = it is the Euclidean 
time coordinate. In fact, all other solutions to the Euclidean equations are divergent [p6| . 
In considering the nucleation of one isolated bubble we may, without loss of generality, 
choose to be real and non-negative. For such configurations the equation of motion 
is 

l_d_ f 3 rf0\ = ldV 

s 3 ds\ S ds) 2 # ' U 

The "bounce" solution [^|] to this equation satisfies the boundary conditions 0'(O) = 0, 
0(oo) = and is characterised by a kink, the bubble wall, localised near s = Ro, where 
Rq is to be determined. For e < 1, i? turns out to be large compared to the width 



: We remark that the choice of potential is not important, and that the phenomena discussed in this 
article are qualitatively similar for any effective potential exhibiting a first-order phase transition. The 
quantitative results will be seen to depend on generic features of the potential, such as barrier height, 
the location of the two minima, and their difference in energy density. 
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of the kink, so one can set s ~ Rq in the region where dip/ds differs significantly from 
zero. Moreover, for small e we have V(4>) ~ Vo(0), where Vo(0) is obtained from V(<f>) 
by setting e = 0, and (f)(0) ~ 77. With these approximations, collectively known as the 
"thin-wall" approximation, Eq. (13]) is easily integrated to yield 



is v W> • (4) 
After one further integration one obtains the bounce solution particular to our potential, 

0b(s) = A= 1 - tanh (r]V\{s - s )) 2 , (5) 

where a constant of integration sq appears because the boundary condition <fi'(0) = is 
not exactly satisfied in the thin-wall approximation. The bounce profile is centred at the 
value of s for which 0b = ^/2, i.e. s = Ro = sq + arccoth(2)/(^ vA) . Near s = Rq the 
bounce solution is well approximated by 

b (s) ~ ^[1- tanh 2fi {s-Ro)} , (6) 



where 

/'o 



-I 0bW = I^) = ^ K ^i (7) 

7] 7] v 8 16 

by virtue of Eq. (Q). Asymptotically, however, the approximation (^) differs from the 
solution (|). Whereas the latter tends to zero as 0b ~ exp(— rjyXs) for s — > +oo and 
tends to r/ as (77 — b ) ~ exp(2?7\/As) for s — > —00, the approximation @ instead ap- 
proaches both these asymptotic values as ~ exp (— 3?7v / A|s|/2). Fig. |l]shows a comparison 
of the analytical approximations of the bounce profile given in Eqs. (D and © with the 
numerical solution of Eq. (|3]). 

The nucleation radius Rq is as yet undetermined, but can be calculated by extremising 
the Euclidean action. For the bounce solution in the absence of gauge fields this action 
can be calculated using the thin-wall approximation and Eq. (f|) without reference to the 
solution itself. One obtains 

S E ~ 2vr 2 (oRl - \ V Ri] , (8) 



r 



where 



a = / ds 
Jo 




d0y / W) = rj 3 VX/2 , V = W ■ (9) 



Here a is the surface tension, or mass per unit area, of the bubble wall, whereas the 
outward pressure p is equal to the difference in energy density between the metastable 
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Figure 1: Bounce profiles for e = 0.2. Here are shown the 
exact solution of Eq. (||) (solid), the thin- wall approximation 
given by Eq. (||) (dotted), and the "linearized" thin- wall ap- 
proximation given by Eq. ([]) (dot-dashed). 



and global minima. The factor 2n 2 R% in Eq. (§) is the area of the three-sphere, whilst 
2ti 2 Rq/4: is the four- volume of the instanton. The action S E is extremised by the value, 

i2 ° = - = T i 7T- ( 10 ) 

V 2er]V\ 

The initial conditions for the evolution of a bubble is given by the t% = slice of 
the instanton, along which d<pb/dtE = 0, such that the bubble is initially stationary. The 
future evolution is found by analytically continuing the bounce solution to Minkowski 
space; details may be found in Ref. f24|. 

The thin-wall limit allows us to make statements about the evolution as well as the 
nucleation of bubbles. A measure of the bubble- wall thickness is provided by fi^ 1 , where 



1 

fi= — 



r) dr 



r(t) _ no -_Mt) (u] 

— ^0-5— — —f. =f , r = — — [LI 
s=Ro Ro VI -r 2 dt 



showing that the thickness of the bubble wall is subject to the expected Lorentz contrac- 
tion. The thickness at t — is given in the thin- wall approximation by Eq. (^). By 
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means of Eqs. ([?[) and fllOD the condition for the thin-wall approximation to be valid can 

be stated 

IM>Ro = — » 1 • (12) 
loe 

2.2 Spherical Membrane Dynamics 

Bubbles whose time evolution is not 0(1,3) symmetric are nucleated, but with lower 
probability, as the result of tunnelling via instantons that do not extremise the Euclidean 
action. In particular, such bubbles may deviate from spatial spherical symmetry, or they 
may be spherically symmetric but not invariant under boosts. We shall here consider the 
latter possibility. 

When the bubble wall is thin compared to the radius of the bubble, it can be approx- 
imated by an infinitely thin membrane. The equation for the radius r(t) of a spherical 



membrane subject to an outward pressure can be written |25] 



dE „ a 9 1 o , , 

^r = °- Bs 7T^7> r ~r • (13) 

where r = dr/dt, a is the surface tension, or mass per unit area, of the membrane and p is 
the difference in energy density between the false vacuum exterior to the bubble and the 
true vacuum inside. Equation (|T3|) simply expresses the conservation of energy, but can 
also be obtained as a special case of the motion of a relativistic membrane in an external 
force field. Because of the relativistic factor (1 — r 2 ) -1 / 2 the first term of E comprises the 
total energy of the bubble wall, mass as well as kinetic energy, while the second term is 
the total energy of the true vacuum inside the bubble. The general solution r(t) is given 
implicitly as 

3E + pr 3 , . 

/ dr = ±t . (14) 

Mo) J(SE + pr 3 ) 2 -9a 2 r 4 



In the special case E — 0, one obtains an explicit solution with 0(1,3) symmetry, 

r 2 - t 2 = Rf , Ri = r(0) . (15) 

If the bubble is nucleated at rest at t = 0, we have from Eq. (|13D that its initial radius 
is Ri = Rq = 3<T I p. Equation ( fT5|) then coincides with the equation of motion for the 
centre of the bubble wall in the 0(1,3) Wick-rotated bounce solution considered above. 

A set of bubble-wall trajectories for different values of E are depicted in Fig. 0. The 
bubbles are all chosen to start from rest at t = 0. If the initial radius Ri is smaller than 
the critical radius i? crit = 2cr/p = (2/3) R , the bubble recollapses, whereas if Ri > R CT it, 
it expands. Only the trajectory with Ri = R is a hyperbola, corresponding to the 
0(1,3) symmetry. It should be noted that, although the initial expansion for i? ; > i? crit 
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Figure 2: Evolution of the radius of spherical bubbles nucle- 
ated at rest in a first-order phase transition. R = 3a /p. 

takes place on different timescales, once the bubble radii have increased by about Ro/2 
the remaining evolution is quite well approximated by a time-translated 0(1,3) solution. 
Therefore, if the average distance between nucleation sites is larger than about 3Rq, the 
0(l,3)-symmetric solution will generally provide a good description of expanding bubbles 
at the time of their collision. 

3 Bubble Collisions and the Kibble- Vilenkin Approx- 
imation 

We now proceed to study collisions of two expanding bubbles of true vacuum. As was 
argued in Sec. |27l], for typical bubbles the phase of the complex field inside each bubble 
is (covariantly) constant, and we may initially choose a gauge in which the gauge vector 
potential vanishes everywhere. Then there will in general be a random initial phase 
difference A9 between the two bubbles. 

It was shown by Kibble and Vilenkin |Tl| that, because of the phase difference, a 
gauge-invariant current jk = ie[(j)*Dk(j) — (Dk(j))*(j)] develops across the surface of inter- 
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section of the two bubbles as they collide. The current gives rise to a ring-like flux of 
the Abelian field strength Fy which takes the shape of a girdle encircling the bubble 
intersection region. Moreover, these authors showed that, under the assumptions that 
the radial mode of the scalar field is strongly damped and a step function can be used to 
approximate the initial phase of the complex field at the moment of collision, the phase 
oscillates and equilibrates to a uniform value as the bubbles merge. 



We shall here review briefly the approach of Ref . |TI[ , as it forms a common basis for 
the analytical methods discussed and developed in this article. We will analyse the results 
arising from their approximation and derive closed-form expressions for gauge-invariant 
quantities such as the Abelian field strength. Crude estimates are extracted for the peak 
value of the field strength and the scale of spatial variation of the field. These quantities 
are relevant for determining the strength and correlation length of magnetic fields that 
may be created through bubble collisions in the early universe in more realistic models. 

Let us consider two non-overlapping bubbles, each locally 0(1,3) symmetric with 
respect to a frame with origin at its centre. The separation of their nucleation events 
is then necessarily spacelike, and one may pick a reference frame and coordinate origin 
such that they are nucleated simultaneously with centres at (t, x, y, z) = (0, 0, 0, ±-R c ). In 
this frame, the bubbles have equal initial radius R\ = Rq. Their first collision occurs at 



(t c , 0, 0, 0) when their radii are R c and t c = y R 2 — Rq by virtue of Eq. (|l"5|). The system 
of two bubbles has the invariance 0(1,2) reducing the problem to a 1 + 1 dimensional 
one that can be expressed in terms of the coordinates z and r, where 

r 2 = t 2 -x 2 -y 2 . (16) 

In these coordinates the point of first contact z = 0, r = t c (denoted by C in Fig. ||b) 
represents the worldsheet of the circle of most recent intersection of the bubble walls in 
3 + 1 dimensional spacetime. This circle recedes superluminarly from the z axis according 
to the equation x 2 +y 2 = p 2 = t 2 —tl, t > t c . Its position at a particular time t is indicated 
by the label C t in Fig. ||a. 

Let us write the scalar Higgs field in polar form = Xe ld . In the axial gauge A z = 
its modulus, phase, and the vector potential take the form 

X = X(t,z), 6(x)=6 a (r,z), A a (x) = x a a(r, z) , (17) 

where a G {0, 1, 2}. The condition 6 l a (r, 0) = fixes the gauge completely. 

Following Ref. |11[ the modulus X of the scalar field inside the bubbles and in their 
overlap region is approximated by the constant true-vacuum value rj. The equations for 



2 Unlike Ref. we do not assume that t c /R c = 1, which would be approximately true only if 
R c S> Rq- 
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-Ft R 

c c 



Figure 3: Bubble collision in two different coordinate systems. 
Dotted curves indicate the extrapolated position of each bub- 
ble wall in the absence of a collision. The union of causal 
futures of bubble- wall intersections, / + (C), extends out to 
the causal horizon I°(C) drawn as a dashed curve, and is the 
only region where physical fields that result from bubble col- 
lision events can exist. Further notation is explained in the 
text, (a) Slice of bubble collision for constant time t and 
constant azimuthal variable (p. Here C t denotes the instanta- 
neous position of the circle of most recent intersection, whose 
radius is p c . The contour ABDEA is used for evaluating the 
Abelian flux, (b) Reduced spacetime diagram of entire bubble 
collision. The point C represents all intersections of the two 
bubble walls. 
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together with 



9g,(r, z) and a(r, z) then reduce to the linear equations 

d% + -d T e & -d% + m 2 e & = o, (is) 

T 

4 

d 2 a + -d T a - d 2 z a + m 2 a = , (19) 

T 

2 
Tfl 

3a + rd T a = — 9 & , (20) 

e 

where m is the vector-boson mass, m 2 = 2e 2 r] 2 . The const ant- modulus approximation 
is not quite consistent, as nontrivial solutions of Eqs. ( |18|) and ( |T9"D do not satisfy the 
remaining equation 

□X-(D M MW)X+~ = (21) 

for constant X = r]. On the contrary, non-trivial gauge-invariant phase gradients 
D^9 = d^9 + eA^ appear inside the causal future of bubble- wall intersections / + (C) 
= {(z, t) I t — t c > \z\) (see Fig. |^b) and serve as sources for fluctuations in the Higgs- 
field modulus. Fortunately, for a small initial phase difference A# the fluctuations in 
X are also small, so that Eqs. flI8|), ( Jl9| ) and ( p0|) remain a good approximation for the 
phase and the vector potential. 

In this article we shall be interested mainly in two quantities, the Abelian field 
strength F az = x a d z a(r, z) and the phase 9 of the complex field. The field strength 
is gauge invariant, but the phase is not. Because the analytical treatment simplifies 
in the axial gauge, while the lattice computations of Sec. |5] are conveniently performed 
in the temporal gauge, it is useful to construct a gauge-invariant generalisation of the 
phase 9 which can be calculated in either gauge, so that the independent results may be 



compared. Such a gauge-invariant phase can be defined [1 1| at any point P by means of 
the line integral of the gauge-invariant phase gradient _D M # = 8^9 + eA^ from the origin 
to P. Unfortunately, this line integral is path dependent. In this work we shall fix the 
definition of the gauge-invariant phase 9 by specifying a path from the origin as follows: 

f(t,x,y,0) r(t,x,y,z) 

6(t, z) = / dx a D a 9 + / dzD z 9 , (22) 

-'(0,0,0,0) J(t,x,y,0) 



where a = 0,1,2 is summed over. By choosing initial conditions for Eqs. fllq ) and ( |19D 
(or the full, nonlinear equations) that adhere to the symmetry of the problem, one can 
arrange for the solution to be such at all times that d a 9 and A a are both odd functions 
of the variable z. Then the first term of Eq. (p2[) does not contribute, and we have 



9(t, z) = [ {T ' Z) dzD z 9(T, z) . (23) 

J(t,0) 
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The expression simplifies further in the axial gauge A z = 0, where one finds 6(t,z) = 
9 a {T, z). Note that, while D^9 is defined locally, the phase 9 is a non-local quantity by 
construction. Therefore, it is not a physical variable, and the evolution of 9 exhibits non- 
causal correlations. Nevertheless, 9 provides a useful, gauge-invariant characterisation of 
the field configuration. 



3.1 Step- Function Approximation 



In the approach by Kibble and Vilenkin [|IT| , a sign function is usedQ to approximate the 



initial phase difference between the two bubbles at the instant when they collide, r — t c . 
This would correspond to bubbles with infinitely thin walls. The initial conditions on the 
phase and its derivative are taken to be 9 a (t c ,z) = 9 sgn(z), corresponding to a phase 
difference of A9q = 2#o, and <9 r # a (t c , z) = 0. The vanishing of the derivative can be 
understood from Fig. |^b, since at r = t c points with \z\ > are in the interior of bubbles 
with constant phase ±# , and 9 a (r,0) = from the gauge condition. In this ansatz, the 
bubbles are treated as though they extend to z — ±oo, whereas the solutions should 
strictly speaking only be trusted where the constant modulus approximation is good, i.e. 



in the interior of the bubbles given by \z\ < R c + y Rq + t 2 . In practise, this does not 
matter, since gauge-invariant physical fields only propagate inside the causal future of 
bubble- wall intersections I + (C). This region extends beyond the overlap region found by 
extrapolating the original bubble trajectories (dotted lines in Fig. |3|), but never exits the 
union of the interior of the two bubbles. 

The vector potential is chosen to be zero initially, a(t c , z) = 0, and its derivative is 
determined from Eq. (|20|). One obtains d T a(t c ,z) = m 2 9 sgn(z)/( y et c ). The solution of 



the linearised equations ([T8|) and ( |19|) for r > t c then becomes |l 1 



. , x 9 t c r°° dk . . ( , , 1 . , . , 
9Jt,z) = / — sin kz costufr — t c ) H sinu t- t c ) ) , (24) 

TIT J-oo k V Ojt c 



. . 9 m 2 t c r°° dk . . 

a(r, z) = — / — sin kz 

Tier 6 J-oo k 



T-t 



— — cosa>(r - t c ) 

+ [- + — t— ) sinu;(r - t c ) 



(25) 



,L0 UJ 3 t c 

where uo 2 = k 2 + m 2 . These integrals are highly oscillatory, but can be put in a form that 



3 A similar approach is adopted in Ref. |lC 
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greatly facilitates their numerical evaluation. Letting T = r — t c we have 



e a (r,z) 



9^ [ Z du ( U(mVf^) - mT MmVT^) 



T JO 



tr 



tr 



VT~ 2 



9q — sgn(z) ( cos mT H sin mT 

mtr 



Similarly for the vector potential A a = x a a we obtain 
m 2 t 



a(r,z) = 6» ~ sgn(z) 

or, for general z, 



cos mT + ( 1 — ! sin mT 



m 2 t 



m m 3 t c 



a(r,z) = - duB^ir.u) 
p Jo 



, \z\ <T-t c 

, \z\ > r - t c . 

(26) 

\z\>r-t c (27) 
(28) 



where p = \fx' 2 + y 2 is the distance from the z axis, and the azimuthal Abelian field 
strength B v = F zp = pd z a has the closed-form expression 



9 m 2 t c 

P 3~~ 

er 6 







rJ§{m\/T 2 — z 2 ) + 



yr 2 - z" 1 

mt r 



J x {m^T 2 -z 2 ) , \z\ < r-t c 



\z\ > T — t r . 



(29) 

The corresponding "electric field" E z is axial and is given by E z = F 0z = —td z a = 
~(t/p)B' p (T, z). The gauge-invariant phase gradient D^9 can also be written in closed 
form; in particular D z 6 = d z Q & is obtained directly from Eq. (p6"l). For a G {0, 1,2} one 
finds, using Eqs. flT|) and 



D n = I 



x a z^r mT —J (mVT 2 - z 2 ) 



VT 2 - z 2 



tr 







\z\ < T — t r 



\z\ > T — t c 



(30) 



The physical fields B v , E z , and D^8 are nonzero only in the causal future of bubble- 
wall intersections / + (C) = {(z, r) \z\ < r — t c }. In the 3+1 dimensional coordinates at 
a particular time t this corresponds to a region enclosed by two rotational paraboloids 
shown as dashed curves in Fig. ^|a. The physical fields are generated by sources on the 
superluminarly expanding circle C t of most recent intersection (p = p c = ^t 2 — t 2 , z — 0) 
and propagate towards the z axis, i.e. into the bubble overlap region. 



15 



One characteristic of the step-function approximation is that the original discontinu- 
ity of 26o in the initial condition for the phase splits after the collision into two equal 
discontinuities of magnitude 6ot c /r at z = ±(t — t c ). Similarly, the Abelian field strength 
has a discontinuity on I°(C); the limit from the interior is Bf z ,,_ t = 6 m 2 t c p / (er 2 ) , 
while it is zero outside / + (C). 

On the plane z = the field strength B v oscillates about zero with an amplitude 
which is an increasing [decreasing] function of p [r]. The largest field strength anywhere 
is found on the circle C t of most recent intersection (p = p c , z = 0) and has the value 
|-Smax| — 0om 2 p c /(et c ). As this circle expands, |-B max | oc t ~ r(t) at late times t, where 
r(t) is the bubble radius. The field strength B v around the maximum defines a peak, 
the radial width of which, Ap, can be estimated from the location of the first zero of 
B^ij, 0) for t > t c (i.e. p < p c ). With m(r — t c ) = Xq at this zero one finds 



= p c - Jt 2 - (t c + — 2 , 31 

V m 2mH 

where xqi < xq < Xu and x n \ is the smallest positive zero of the Bessel function J n (x), 
xoi ~ 2.405, xn ~ 3.832. The last expression in Eq. (^l|) is the leading behaviour at 
large times t. The narrowing of the peak with time is a relativistic effect. For a fixed 
value of r, the field-strength peak extends in the z direction across J + (C), so that the 
width in that direction can be estimated as Az = Xq/tti. 

We can now derive a crude estimate of the Abelian flux contained in the maximal 
peak by taking 

#0 2/1 , x o 



$B, P cak~ |S max |ApA^~^x 2 ^l + ^) , (32) 



where the last expression is the constant flux approached at late times. One can compare 
this to the total flux of B v through a semi-infinite plane based on the z axis (i.e. the 
total flux of all field lines that "run around the z axis"). This flux, $s, is obtained by 
evaluating the line integral of at constant time t > t c along the path ABDEA indicated 
in Fig. |3]a. In the axial gauge A z = 0, the segments AB and DE do not contribute to 
the integral. In order to enclose all flux of the Abelian field strength, the segments BD 
and EA must be chosen so that their z coordinates satisfy \z\ > t — t c , and the points 
D and E must be located at p > p c (i.e. r < t c ). Assuming that the vector potential A k 
vanishes for r < t c , in accordance with the initial condition, only p such that < p < p c 
contribute to the integral. From Eq. (|27|) we find 



= -([ +[ )dx a A a = 2[ drra(r,\z\) 

, t>t c . (33) 



29o 

e 



1 ( cos m(t — t c )) H sinm(t — t c 

t V rate 
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Note that the integral receives no contributions from regions where the Higgs-field mod- 
ulus X is fluctuating, and therefore the expression ( p3[) for the total flux is valid also for 
large initial phase difference A9q of the bubbles (see discussion after Eq. (|2ID). In this 
sense, the total flux is of a topological nature. 

As a function of time, the Abelian flux rises from zero and begins to perform 
damped oscillations about the late-time limit A#o/e, where A6 = 2# is the initial 
phase difference of the two bubbles. The amplitude of these oscillations decreases as 
We emphasise that the flux does not go away; it is there to stay. 

Comparing with Eq. (j32l) we see that the flux <&s,peak contained in the maximal peak 
is of the same order of magnitude as (or somewhat larger thanF]) the total flux This 
shows that a significant portion of the total flux is contained in an outward travelling 
circular flux tube near p = p c , z = whose field strength -B max is increasing with time. 



4 A New, Smooth Analytical Solution 

The step-function approximation described in the previous section gives a rough picture 
of the qualitative behaviour of the physical fields in a U(l) bubble collision. For instance, 
by looking ahead to Fig. |] one can see that the shape of the profile of the Abelian field 
strength in the bubble overlap region is quite well reproduced. 

On the other hand, the step-function approximation predicts a peak value of the field 
strength that differs by a factor of four from the correct value obtained in numerical 
simulations (a factor of sixteen difference in energy density). Moreover, the solutions 
(p6l)-(l30D have discontinuities on the future null surface of bubble intersections I°(C) 



and are thus unable to account for the behaviour of the fields near that surface. Finally, 
the solutions are limited to the region r > t c (where r = y/t' 2 — p' 2 ) and therefore provide 
no clue to the field behaviour shortly before the bubble collisions or, at later times, 
outside the bubble walls. 

Realistic bubble walls of finite width, such as given by Eq. (||), have exponential tails 
that begin to overlap long before the instant of collision. The overlapping tails initiate the 
evolution of the phase of the scalar field and of the vector potential so that, in contrast 



with the assumption made in Sec. |3.1] , the functions d T 9 a and A a = x a a are already 
nonzero when the collision reaches full impact at r = t c . 

In this section, we shall model the initial evolution of 6 a and a and obtain a smooth 
analytical solution. To do so, we first observe that the scalar Higgs field of two well- 
separated bubbles can be written [II], [L3| 



= $ x (x) + $ 2 (x) , (34) 



^This is not a contradiction, because B v changes sign in the bubble overlap region 
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where $1, $2 are 0(l,3)-symmetric solutions of the Higgs field equation with $j — > 
away from the bubbles and approaching the vacuum expectation value 77 well inside 
each bubble. In the absence of gauge fields the linear superposition remains a good 
approximation even when the bubbles begin to overlap, because in the region of overlap 
the Higgs field is small and the field equation can be linearised about $ = 0. When 
gauge fields are included, the reliability of the linear superposition approximation for 
describing the phase of the Higgs field depends on the choice of gauge. We have found 
that in the temporal gauge the initial evolution of the phase is well described by a linear 
superposition approximation, as long as the bubble overlap is weak. 

The position of a bubble wall can be defined as the set of points where = rj/2, i — 
1,2. Note that this definition gives the location of each bubble wall in the absence of 
the other bubble. The bubble collision is then defined as an event x for which |$i(x)| = 
1^2(^)1 — 7 ?/2. This corresponds to the intersection of the extrapolated trajectories of 
two well-separated bubble walls. When the bubbles begin to overlap strongly, the actual 
scalar field <&(x) has a nonlinear behaviour and therefore cannot be used to define the 
event of collision. 

Using the same geometry as in Sec. 0, depicted in Fig. |3|, the assumptions underlying 
the present analytical treatment can be formulated as follows: The bubble collision occurs 
at t = t c , z = with r = y/t 2 — p 2 . For r < t c the field evolution is given by the linear 



superposition (p4|). This defines the fields outside the bubble walls as well as before the 
time of collision. For r > t c the scalar field is assumed to have constant modulus |$| = rj. 
This reduces the problem to one of finding the solution of the linearised equations ( |TBD 
and fll9l ) with initial conditions at r = t c provided by the linear superposition. 

Note that although linear superposition is an excellent approximation for r < t c , 
and the const ant- modulus approximation works very well for r ^> t c , both approxima- 
tions break down near r = t c . The approach taken here is therefore one of joining two 
asymptotic solutions across a non-linear region where no solutions are known. 

Some further approximations are necessary. We assume that the thin-wall approxi- 
mation (1T2| ) is valid. Moreover, we shall require that the two bubbles overlap negligibly 
at the time of their nucleation, t = 0. At this time, the bubble walls are a distance 
2(R C — Rq) apart, and the thickness of the bubble wall is given by Pq 1 (refer to Eqs. (^) 
and (|TT|)). One obtains the no-overlap condition 

2p (R c - Rq) > 1 . (35) 

Because of the exponential decay of the bubble-wall profile, results accurate to about 
3% are obtained already when the left-hand side of Eq. (|35| ) is greater than 3.5. Via the 
identity t 2 = R 2 — R% the condition of no overlap also implies the following very useful 
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relation, 



2(3 = 



Rc 



< 1 . 



(36) 



2^t\ 

Let us now consider two bubbles, one with constant Higgs phase 9 centred at x — 
y = 0, z = R c and the other with constant Higgs phase — 9 Q centred at x — y — 0, 
z = —R c - It is convenient to express the Higgs modulus in the bubbles as a product of 
the vacuum expectation value r\ and a function / of the variables 

[z T Rc) 2 + r 2 



Rq 



2R r 



(37) 



which possess 0(1,3) symmetry with respect to the nucleation events (0, 0, 0, ±R C ), re- 
spectively. The modulus of the Higgs field in each bubble is then 



vfM = M^ 2 + v 2 + (zt Rc) 2 - 1 2 



Rq — 2R c u± 



(38) 



where 0b is the approximate bounce solution given by Eq. ([|) or, alternatively, Eq. (Q). 
The positions of the bubble walls correspond to u± = 0. The profile function f{u±) has 
non-negligible variation only near the bubble walls or, more precisely, for \u±\ < 
where /z" 1 = Rq/(hqR c ) is the width of the bubble wall at the collision time t c given by 
equation ([11]). Using the thin- wall condition ( |T2| ) it is easy to show that /i" 1 -C Rq/{2R c ). 
The argument of the bounce profile, Eq. ([38]), can therefore be safely linearised, giving 

(39) 



rif(u±) w b (i? o - R c u±/R ) . 
The linear superposition (|34|) can now be written 



$(x) = X(x)e 



i0 t (x) 



V 



f(u + )e Wo + /(«_)e" 



We obtain the following expressions, which are assumed to be valid for r <t c : 



X 2 
Ot 

o T e t 
d z e t 



V 2 {[/(«+) + f(u-)] 2 cos 2 9 + [f(u+] - f(u-)] 2 sin 2 £o} , 
/(«+) - /(«_ " 



arctan 



sin 29 f] 



tan 9p 



/(«+) + /(«_) 



2 



sin 2# ^ 



[/'K)/(r)-/'( U -)/W] , 



7? r 



(40) 

(41) 
(42) 

(43) 
(44) 



where 6 t is the phase of the Higgs field in the temporal gauge. 
At the collision the argument u± simplifies to 



u± 



±z it 



2R r 



T = t c 



(45) 
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We see that u± is zero not only in the plane of collision z = but also on the two remote 
bubble walls located at z = ±2i? c . However, no interesting dynamics takes place near the 
remote bubble walls, and we can ignore their presence. This is because physical fields 
excited by the collision can extend at most a distance ~ /i" 1 outside the null surface 
7°(C) drawn in Fig. |3|. The protrusion beyond the ideal null surface is due to the finite 
thickness of the bubble wall. It can be shown that the distance between I°(C) and the 
remote wall remains always larger than R c + t c . Therefore, because of the no-overlap 
condition fl35f) and the thin- wall condition (|i~2"D , physically interesting fields will never 
reach the remote bubble walls. 

Concentrating instead on the collision region near z = 0, we notice that, for r = t c , 
the profile function f(u±) has non-negligible variation only for \z\ < /i" 1 and quickly 
approaches a constant value for larger \z\. From Eq. fl35|) we find that p^ 1 <C 2R C and 



we can thus ignore the term of order z/R c in Eq. (45), obtaining 



u± « ±z , t = t c . (46) 

Similarly, for r <t c the products f'(u±)f(u T ) in Eq. flip) are non-zero only for \z\ < /x" 1 , 
and we can neglect the term z/R c also in this equation, which reduces to 

2 

d z 9 t = sin 29o ^ [/'(«+)/(«_) + /'(«_)/(«+)] , r < t c . (47) 
4.1 Solution for the Abelian Field Strength 

In this section, we show that the initial evolution of the scalar Higgs field in the linear 
superposition approximation generates non-zero vector fields already before the collision 
time t c from the exponential tails preceeding the expanding bubble walls. The vector 
field is evolved to r = t c , and its values there are used as initial conditions for the ensuing 
linear evolution. 

For non-constant Higgs modulus X, the z component of the Maxwell field equation, 
Eq. (PU|), generalises to 

d z (3a + rd T a) = 2eX 2 d z 6 SL , (48) 

where 9 a is the phase of the Higgs field in the axial gauge defined by A 3 = 0. This phase 
can be expressed in terms of 9 t , the phase in the temporal gauge defined by A = 0, as 
follows 

2 

a = 0t-| T Cid(T 2 ) , (49) 

where we used the initial condition that N 1 — > and therefore # a — 9 t —>■ as r 2 = 
t 2 — p 2 — ► — oo. One obtains the differential- integral equation 



d z (3a + rd T a) = 2eX 2 (d z 6 t ~ ~ f d z ad(T 2 )j . 



(50) 
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We shall show at the end of this section that the integral term on the right-hand side of 
this equation is subdominant and can be neglected provided that 

f3 2 m 2 t 2 c cos 2 6 ^ < 1 . (51) 
Using Eq. (K7I) the equation for a then becomes 

2 

77? 

d T (r 3 d z a) = — sin2fl r9 T [/(«+)/(«_)] . (52) 
e 

We now specialise to the approximation of the bounce solution given by Eq. @. 
Using Eq. ( |39D one obtains 

/H = ^[l + tanh(2/i c n)] . (53) 

Equation ([52]) can then be integrated provided that we make another well founded approx- 
imation. Realising that the composite function (f ou)(r) is monotonically increasing for 
all z, we have for r <t c the leading behaviour f(u) ~ g(z) exp[4// c i c (r — i c )/.R c ] for some 
function g(z). The right-hand side of Eq. (0) consequently has the leading behaviour 
~ h(z)r 2 exp[8/j c t c (r — t c )/i? c ] for some function h(z). Its integral with respect to r is 
thus dominated by contributions from an interval t c — 6 < r < t c , where 5 = R c /(8fi c t c ) = 
(3t c /2 t c . Let us now consider the identity r 2 = t 2 c + 2t c (r — i c )[l + (r — t c )/(2t c )}. 
Near r = t c we can replace r 2 by its linear expansion, provided that \t — t c \/(2t c ) 1. 
However, the latter condition is always satisfied in the interval (t c — 5, t c ] by virtue of 
Eq. (p6|). We can therefore consistently replace the variables u± by new variables v±, 
linear in r, defined by 

v ± = ±z + ^(r-t c ) . (54) 

A term -z 2 /(2R c ) was omitted from v± in analogy with the discussion leading to Eq. (|4q) . 
To justify this, let us again ignore the remote bubble walls located at z = ±(R C + 
J 'Re + t 2 — t 2 ) and take into account that f{v±) has non- negligible variation only for 
\v±\ < /j" 1 . It then follows from the triangle inequality and Eq. fl35"D that \z\ < \v±\ + 
t c \r - t c \/R c < 9fi~ 1 /8 < 2i? c in the interval t c - 5 < r < t c . 

The linearisation ( |54] ) has the added benefit that f(v±) has an exponential fall-off as 
r — > — oo (in contrast, f{u±) approaches unity as r — ► — oo and falls off exponentially only 
for t 2 — > — oo). Because the contribution from r t c — S is negligible, for the purpose 
of integration we can extend the domain of definition of the fields to t 6 (— oo,t c ], 
identifying only positive r with the physical quantity \ft 2 — p 2 . 
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By means of integration by parts, and using the initial condition a(r 
the solution of Eq. fl5"2]) is written 



-oo, z) 



d z a 



m 2 R c 
e t 3 
m 2 R r 



sm29 \rf(v + )f(v. 



f{v + )f{v_)dr 



sin 2#o{t(i + tanh 2/i c t> + )(l + tanh 2/i c 



e 4r 3 

- 2(3t c [ln(l + e 4/ ^+) + ln(l + e 4 ^" 
+ 2(3t c [ln(l + e 4 ^+) - ln(l + e 4 ^" 



coth4/i c 2;| 



(55) 



where (3 is defined in Eq. (0). The expression for d z a can be integrated, using the 
condition a(r, 0) = prescribed by symmetry. One obtains 



a{j, z 



sm20„ 



rfe 2 



r 3 e 



In 



1 + 6c 



6 2 - 1 b + c 



ln(c)ln(l-6 2 ) + Li 2 
b(l + c)\ 



Li 



1-b 



Li 2 



6(1 + c) 
c(l - b) 
' b{l-c) \ 
l + b J 



I + Li 2 
+ 2Li 2 ( 



c(l + b) t 
-be) 



2Li 2 



where 



exp 



tic 



c = exp(4/i c z) 



and where the poly- logarithm function Li n (z) is defined by 

(-1)" 
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U n (z) 



[n 



dx 



In 71 " 1 x 



oo z k 



X 



l/z 



E 

k=l 



k' 



(56) 



(57) 



(5? 



For real z and n = 2 the power series is convergent for |z| < 1. The evaluation of Li 2 (z) 
for z near or outside the convergence radius of the series is discussed in Appendix B. 

The expressions (^) and ( pB[ ) are valid for t c — 5 < r < t c and approach zero rapidly 
for r < t c — 5. Nevertheless, these expressions for d z a and a have poles at r = 0. The 
poles are an artifact of the linearisation (|54D and would cancel in the exact solutions, 
leading to finite (and exponentially small) functional values of a and d z a as r — > 0. 
Because 5 <C t c one can easily remedy this slight problem by explicitly setting a(r, z) = 
for r <C t c — 5. 

As the next step, we shall now find expressions for a(t c , z) and d T a(t c , z) to use as 
initial conditions for the linear equations governing the evolution for r > t c . Taking the 
limit r — > t c in Eq. (|56"D, using the asymptotic result ( B.5|) for the dilogarithm Li 2 (a;) 
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(see Appendix B) as well as the relations 
m 2 sin 29 1 



a(t c , z) = (3 



and ( |B.3| ), we find 
tanh 2fi c z + (3 8/x c zln(2 cosh2yU c z) 



e 2 

+ Li 2 (tanh2 / u c z) - Li 2 (- tanh 2{jt c z) + 2 Li 2 (-e 4 ^ 2 ) - 2 U 2 (-e~^ cZ ) 1 1 (59) 



7r 2 \ m 2 sin 2#n 



tanh 2fi c z 



(60) 



The expression ( |60| ) is an excellent approximation; it was obtained by matching the 
functions in the limit of large \z\ using Eq. ( |B.5| ) and the special values Li 2 (0) = 0, 
Li 2 (-1) = -tt 2 /12 and Li 2 (l) = vr 2 /6. 

Similarly, by differentiating Eq. (|56| ) with respect to r and taking the limit r — ► t c , 
we obtain 



d T a(t c , z) 



3 , . m 2 sin 29 

—a(t c , z) + — — tanh2/i c 2; 



7T 



1 - 3/5 + -f3 2 



m 2 sin 2# 
etT 2 



tanh2/z c z . 



(61) 



Note that, in the limit of an infinitely steep bubble wall (/i c — > oo, /5 — ^ 0) and for small 
Higgs phase difference 29o, Eqs. (|59|) and (plf) reduce to the initial conditions a(t c , z) = 0, 
d T a(t c ,z) = m 2 9oSgn(z)/(et c ) of the Kibble- Vilenkin step function approximation (see 
Sec. |pand Ref. 0). 

The solution of Eq. ( |I9|) for r > t c subject to the initial conditions (|60|) and ([H]) is 
expressed as 



a(r, z) 



m 2 sin26 l ot c Z" 00 sin&z 



5A*c ^ 



sinh ^ 



-^+/3(l-^)t c r 
a^t,. 12 



cosu;(r — £ r 



+ 



,7T 2 ,i, 



CO tuHr 



sincufr — £ r 



(62) 



where uo 2 = m 2 + k 2 . We have retained as many orders of the small parameter {3 as is 
necessary to ensure that this solution joins smoothly with the solution for r <t c given 
in Eq. (|56|) . 

The Abelian field strength is defined by B v = pd z a, where p = yfx 2 + y 2 and where 
d z a(r, z) is given by Eq. (]55f) for r < t c , or obtained for r > £ c by differentiating Eq. (|62]) . 



4.2 Solution for the Gauge-Invariant Phase 

Let us now turn to the gauge-invariant phase 9. Because of Eq. (|23|) it is identical to the 
phase 9 a in the axial gauge defined by A 3 = 0. The phase # a is given by Eq. (pE9|) which 
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has two terms, one being the phase 6 t in the temporal gauge and the other an integral 
of the vector-potential field a(r,z). For r < t c , these can be evaluated using the linear 
superposition approximation. 

The phase 8 t is given explicitly by Eq. ( f4"2"D . For r < t c we shall therefore use the more 
complicated, but precise, expressions for the bubble profile function with f(u) given by 
Eq. ( |3"%| ) and </>b given by Eq. (^). This bounce solution b has the correct exponential 
tail far from the bubble wall, which is particularly important for the global behaviour of 
gauge-invariant phase 6. 

The integral in Eq. (f49|) can be evaluated using the expression for a(r, z) given in 



Eq. (pop , provided that one expands inverse powers of r in a Taylor series about r = t c . 
It suffices to do the calculation to lowest order in the parameter (3. One finds 



dr ra = j3 2 m 2 t : 
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(63) 



which is valid for r <t c . In particular, at r = t c the expression reduces to 

el drra = (3 2 m 2 tl — — — |4/x c z ln(2 cosh 2/i c z) + Li 2 (tanh 2/i c z) — Li 2 (— tanh2/i c z) 
J—oo 2 L 



+ Li 2 (l - e 4 ^ 2 ) - Li 2 (l - e~ 4 ^ 2 )) - Li 2 



1 + e 4 ^ z 



+ Li 2 



1 



1 + e ~ A ^ z 



f3 2 m 2 t 2 



sin 29 7r 2 
2 12 



tanh2/i c 2; . 



(64) 



We are now ready to derive approximative expressions for a (i c , z) and d T 9 a (t c , z) to 
be used as initial conditions for the ensuing linear evolution. First, from Eqs. fl4"2] ) and 
([46|) we obtain 



Bt(t c ,z) 



arctan 



tan 8q 



m - f{-z) 



'f(z) + f(-z) 
arctan(tan6 , otanh2/i c z) « ^otanh 



tan6*n 



0n 



-2// c 2 



(65) 



where we used Eq. (|53|) in the intermediate step. Substituting this result together with 
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Eq. 



into Eq. 



we find the initial condition 



0&(tc,z) » 6> tanh 



tan 6*0 



-2/i c z 



(3 2 m 2 t 2 



sin 26 Q tt 2 



12 



tanh2/i c z . 



(66) 



Turning next to the r derivative given by Eq. ([431), using Eqs. 

2fi c t c 



R r 



tasx 9 g(fi c z) 



j) and (|53D we have that 
(67) 

(68) 

cosh 2x + tan 8q sinh 2x 
We shall need an approximation of the function g which possesses a simple Fourier 
transform. To this effect, let us define the value xq by 



d T e t (t c ,z) 

where the function g(x) is defined by 



tanh 2x 



satisfying 



Xq — - COsh 1 



V5 + 4 cos 26»o 



~(1 + v /5 + 4cos2^ ) 



i 

l\ 2 



cos 2 9 



VV5 + 4cos2# + 3/ V5 + 4 cos 200 + 1 + 2 cos 26 
An excellent approximation to g(x) is then given by 

g{x) « 16^g(x ) + . 

Combining Eqs. fl6"T|), ([7T]), ( pCf ) and ( |4~9l ) one finds the other initial condition 

d T 6 a (t c , z) = d T 8 t (t c , z) - et c a(t c , z) 



(69) 



(70) 



(71) 



32fi c t c 



tan 6n 



-P [1-/3 



7T 

12 



m t r 



sin 200 



tanh2/z c z . (72) 



(S* 2 , + £z*) 2 

The solution of Eq. (|r8"|) for r > t c subject to the initial conditions (^) and (|72|) is given 
by 



t f°° „ sinkz 
9Jr,z) = — I dk— x 



r J-oo 
01 



4/i c 



tan # n sinh 



6q 7rk 
tan So 4/i c 

^o 2 



sin 29 f3 2 7t 2 m 2 t 2 c 

2 12sinh (S) 
sin 2# /?m 2 t; 
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2+2 
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tan 



sinh ( 



6>p 7rfc 
tan $o 4/i c 

32t? 



sinh (ft) 



tan flo ^-^ gf(x )a:oA;exp(-v / 3xo|fc|//i c 



sin a;(r — £ c ) 



LOt c 



(73) 
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where u 2 = m? + k 2 . 

Finally, let us justify the approximation made in Eq. (|50"D, where the integral term 
was neglected. To estimate the relative significance of the two terms in the right-hand 
side of that equation, recall that the integral / d,Trd z a is an increasing function of r which 
is significantly nonzero only in a narrow interval t c — 5 < r < t c , where 5 = R c / (8fi c t c ) = 
(3t c /2 -C t c . The largest values of the integral occur at r = t c , where it is easily compared 
to the other term. From Eq. ( fMl) one obtains 

/tc 77"^ 
drrd z a m /3 2 m 2 t 2 c /j, c sin 26 — sech 2 2/i c ;z , (74) 
-oo 12 

and from Eq. Q65Q , similarly, 

d-A{t c , z) « 2/2 c tan 6 sech 2 ^!!±2fi c zj . (75) 
Therefore, the condition for the integral in Eq. ([50]) to be negligible is given by Eq. (|51|). 



5 Lattice Computation 

In order to compare our approximate analytical solutions to the true evolution of fields 
in the original field theory, we have evolved the system of two bubbles numerically on the 
lattice, using the temporal gauge A = and the techniques described in Appendix A. As 
the first step of such a lattice simulation it is necessary to specify the initial conditions. 
In particular, one must determine the initial shape of the bubble walls. 

For definiteness we consider the Higgs potential (p]) with the fields and coordinates 
rescaled so that A = r\ — 1. We choose e = 0.2, corresponding to the thin- wall condition 
(0) being marginally satisfied. This means that the bubble wall is thin, but not extremely 
thin, compared to the nucleation radius Ro (see Fig. 0). The gauge coupling of the 
rescaled fields is taken to be e = 0.5. 

The Euclidean field equation relevant to the nucleation of a single spherical bubble 
at rest, Eq. (|j), has been solved using a shooting method. At t — in Minkowski 
space the spatial dependence of the single-bubble configuration coincides with that of 
the Euclidean instanton solution. The initial two-bubble configuration is approximated 
by the linear superposition of the scalar fields of two well separated single bubbles. The 
vector potential A 1 and the time derivatives of all fields are set to zero initially. 

We have used a spatial Cartesian grid consisting of 125 3 unit cells with a lattice 
spacing of a = 0.2 (this a need not be confused with the field a(r, z) defined in Eq. (|lTD). 
In order to maximise the dynamical range, we placed the bubbles in two adjacent corners 
of our box and employed reflective boundary conditions. More precisely, on each face 
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of the box we imposed dQ/dx 1 - = 0, A 1 - = 0, and dA^/dx 1 - = 0, where _L signifies 
the outward normal and A" denotes the two vector-potential components parallel to the 
face. Derivatives were approximated by symmetric differences, which caused the planes 
of symmetry in each direction to be located at the lattice points with indices 2 and 124. 
Thus, the distance between bubble centres is 2R C = 122a = 24.4. 

Finally, the system was evolved using a time step of ha with h = 0.2. The lattice 
Gauss constraint, Eq. ( A. 18 ), was monitored and found to remain zero at the level of 
machine accuracy, i.e., one part in 10 15 . 



6 Analysis 

We shall here present results for the Abelian field strength and the gauge-invariant 
Higgs phase 6 = J dx^D^O in the collision of two expanding bubbles. Analytical results, 
within the Kibble- Vilenkin step-function approximation as well as our new approach, 
are compared to the true field evolution obtained from the numerical simulation on the 
lattice. 

First, let us review the values of various parameters that occur in the solutions. As 
was mentioned in the last section, we take A = r] — 1 and e = 0.2 as parameters of 
the Higgs potential. The gauge coupling is e = 0.5 and the bubble radius at the time 
of collision is R c = 12.2 . We then find directly the vector-boson mass m = 
0.7071 . Expressions (0) and ( |TUD determine, within the thin- wall approximation, the 
initial bubble-wall slope /iq = 0.375 and the nucleation radius Rq = 7.5 . The time of 



collision is obtained from Eq. (|l~5| ) with R± = Rq, which gives t c = J R"^ — Rq = 9.62 . The 
bubble-wall slope at the time of collision fi c = fi R c /R = 0.61 is found from Eq. (|TT|) . 
As a consequence, (3 = -R c /(4/x c t^) = 0.054 < 1 as required. Finally, the constant so 
occurring in Eq. (|5]) is given by So = 6.95. This completes a list of parameters that are 
uniquely determined by the physical model and by the chosen distance between nucleation 
centres, 2R C . The average value of this bubble separation can also be calculated from 
the model parameters [[Jl[ . 



The remaining parameter, #o> is determined by the random nature of spontaneous 
symmetry breaking. We choose 6 = n/8, corresponding to a particular case in which two 
colliding bubbles have a Higgs phase difference of 7r/4. This phase difference is sufficiently 
small that important approximations remain valid, such as the assumed constancy of the 
Higgs-field modulus inside the bubble overlap region after collision. Moreover, it can 
be checked that the no-overlap condition ( |35D is satisfied with the left-hand side being 
approximately 3.5 ; this ensures that the overlap of two bubbles is negligible at the time 
of nucleation t = 0. Finally, the condition (0), used in the derivation of the analytical 
solution for 9, is satisfied with the left-hand side being approximately 0.095 <^ 1 . With 
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the chosen value of 9 , Eqs. (F9"|) and ( |70| ) give xo = 0.314 and g(a;o) = 0.365. 

The results for the Abelian field strength are most easily explained with the help of 
Fig. The first and third columns of this figure depict results from the numerical evo- 
lution of the fields on a lattice, which was described in Sec. |5|. Because of the cylindrical 
symmetry, the fields depend only on the coordinate z along the symmetry axis and the 
distance p from this axis. Two bubbles were initially nucleated in the far left and near 
corners, defined by p = 0, z = ±12.2, with nucleation radius Rq = 7.5. In column (a) we 
can follow the evolution of the modulus of the scalar field as the bubbles expand. The 
bubbles first collide at time t c = 9.62, which falls between the first and second frames. 
As they continue to expand, the modulus within the region of overlap of the two bubbles 
begins to fluctuate. The fluctuations are moderate, and therefore the approximation of 
the Higgs modulus by a constant value |<&| = r\ within the union of interior of the two 
bubbles is well justified. A larger distance 2R C between nucleation centres would lead to 



higher bubble- wall speed at collision, v c = J 1 — Rq/R 2 , and therefore to larger fluctu- 
ations in the Higgs modulus. Violent fluctuations in the modulus are induced for large 
phase differences A9 > 7r/2 regardless of the collision speed. Therefore, we don't expect 
our analytical treatment to remain valid for large values of 8q. 

Turning next to the field strength B v depicted in column (c), we notice that the 
flux of B v is mainly concentrated in a peak which moves away from the z axis in the 
symmetry plane z = 0. Because of the cylindrical symmetry, the peak corresponds to a 
circular flux tube enclosing the region of bubble intersection. In fact, by comparing with 
column (a), we see that the peak is located near the value of p at which the bubbles have 
most recently come into contact. This circle of most recent intersection is represented by 



the point C t in Fig. ||a and has radius p c = yr 2 — t 2 . As a consequence, the flux tube 
moves away from the symmetry axis with a speed greater than the speed of light. This 
is not in contradiction with special relativity; the field strength in the peak has been 
generated locally and causally by the most recent bubble-wall intersection, while the tail 
toward the bubble interior consists of fields generated by earlier intersections occurring 
at smaller radius p. 

In column (b) we observe that the solution obtained by Kibble and Vilenkin us- 
ing a step-function approximation (see Sec. [Tl] and Ref. [PUD gives a very reasonable 
qualitative description of the behaviour of the field strength in the region I + (C) t = 
{(z,p) | (t c + \z\) 2 + p 2 < t 2 }, which constitutes the causal future of bubble intersections 
and is shown in Fig. ^a. On the other hand, the peak value of the field strength in this 
approximation is about four times higher than the true peak value obtained from the 
lattice evolution. The reason for the higher field strength is that the step-function initial 
condition leads to a much higher time derivative of the field strength at the collision than 
does the smooth initial condition (FID obtained with finite-width bubble walls. 
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Figure 4: Bubble collision with initial Higgs- field phase difference A9 = 29 = tx/A. 
a) Lattice evolution of the modulus |$| of the Higgs field, b) Kibble-Vilenkin approximate 
solution for the azimuthal field strength B v , c) Lattice evolution of B^, d) Our new 
analytical solution for B^. Here p = (x 2 + y 2 ) 1 ^ 2 is the distance from the symmetry axis 
of the collision, |<E>| is measured in units of 77, B v in units of m 2 /e, and p, z, t in units 
of A -1 / 2 ?? -1 . The centres of the two bubbles are at p = 0, z = ±12.2. The bubbles first 
collide at time t c = 9.62, which falls between the first and second frames in each column. 



The new, smooth analytical solution is displayed in column (d). It is found to be 
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in excellent agreement with the lattice results (c).Q The two solutions are more easily 
compared in Fig. [|, which shows the field strength on the symmetry plane z = 0. Here 
the agreement of the new solution with the numerical results is rather remarkable. The 
tail that leads down from the peak towards higher values of p is accurately reproduced. 
We find that the maximum field strength increases with time, and the width of the peak 
Ap decreases with time as predicted by the analysis in Sec. pTT| . Note that there is a 
trough with negative values of the field strength trailing the primary peak by about a 
peak width. This trough corresponds to an expanding circular flux-tube loop with smaller 
radius that the primary one and with the Abelian field strength pointing in the opposite 
direction. Such a secondary flux tube is created by wave-like excitations travelling inward 
from previous locations of the primary peak. In fact, the two-dimensional nature of the 
bubble-collision problem implied by the 0(1,2) symmetry explains why these solutions 
are similar to waves on a water surface. The trough is less accurately reproduced by the 
analytical solution as time progresses. 

Let us now consider the gauge-invariant Higgs phase 9 = J dx^D^O, properly defined 
in Eq. (|23|). This phase is equal to the Higgs phase 9 a in the axial gauge. Fig. |6] shows a 
comparison between the numerical solution for the phase and the analytical solution. As 
the bubbles collide, the phase in the interior of the bubbles begins to oscillate (to see the 
location of the bubbles at each time step, refer to Fig. ^a). Note that 9 changes instantly 
also for large values of \z\, which may at first seem unphysical. The phase 9 is, however, 
a non-local quantity by construction, obtained by integrating the local gauge-invariant 
phase gradient D^9. One can check that D^9 is non-zero only within the causal future 
of bubble-wall intersections, as should be expected from a physical quantity. Although 
9 has no direct physical meaning, it provides a useful, gauge-invariant characterisation 
of the field configuration, by means of which different solutions can be easily compared. 
It is also a very sensitive measure, because small deviations in D^9 accumulate as they 
are integrated. The agreement between the numerical and analytical solutions in Fig. |6] 
is nevertheless quite good. Through careful study one finds that the analytical solution 
is slightly delayed in time after the collision, and also that the amplitude of oscillation is 
smaller in the analytical solution than in the numerical evolution. We believe that these 
small deviations are the result of non-linear evolution near the time of collision that our 
analytical approach is unable to account for. 

A particularly interesting quantity is the gauge-invariant phase difference of the two 
bubbles. For r > t c this is obtained by considering the z — > ±oo limits of Eq. (|73|) , while 
Eqs. (fi2"|), fl63|) and ( fj9|) in the same limits give the result for r <t c . With the definition 

5 In fact, the agreement is even better than it appears to be because, in order to produce black-and- 
white surface graphs such as these, one must use a sparse mesh which cannot always be brought to 
coincide precisely with peaks and other features. 
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Figure 5: Field strength (in units of m 2 /e) on the plane of symmetry z — 0. The 
numerical solution (solid curve) is compared with the Kibble- Vilenkin approximate so- 
lution (dashed) and with our new analytical solution (dot-dashed). For t > 12, the 
Kibble- Vilenkin solution extends to higher values of B^ than what is visible in the graph. 
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a) b) 

Gauge-invariant phase Gauge-invariant phase, 

from lattice evolution new analytical solution 




Figure 6: Gauge-invariant phase 6 = J dx^D^O in a bubble collision with initial Higgs- 
field phase difference A6 = 26 q = it/ 4: a) Lattice evolution of 6, b) Our new analytical 
solution for 6. Again, p = (x 2 + y 2 ) 1 ^ 2 is the distance from the symmetry axis of the 
collision and p, z, t are measured in units of A -1 / 2 ?? -1 . The bubbles first collide at time 
t c = 9.62, which falls between the first and second frames in each column. 
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A(9(r) = 6(t,z -> +00) 
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(76) 



(77) 



where b is denned in Eq. ( p7f) - In particular, on the axis through the bubble centres 
(p = 0) the phase difference is given by the above expressions with the substitutions 
r — > t and b — ► 6 = exp[4/i c t c (t — t c )/i? c ]. 

The oscillatory behaviour of the phase difference A6(t) after the collision (t > t c ) is 
fully accounted for in Eq. ( [76]) . The amplitude of oscillations is proportional to t" 1 . The 
phase difference between the interior of the bubbles therefore equilibrates to zero as time 
progresses. This behaviour was first noted by Kibble and Vilenkin [JTH| . In fact, in the 
limit of infinitely thin bubble walls (/3 —> 0), Eq. ([76{) agrees with Eq. (p6|) , which is their 
result. On the other hand, the phase difference far from the symmetry axis remains 28q . 
This is evident in Fig. |6]a, but also follows analytically from Eq. ([FT]) in the 6^0 limit, 
using Eq. ( |B.4|) of Appendix B. 

Finally, we shall derive the total flux of the Abelian field strength. This is the flux 
through a semi-infinite plane based on the z axis, and can be computed in analogy with 
the discussion leading to Eq. ([!3|) by evaluating the line integral of the vector potential 

at a constant time t along the closed path ABDEA indicated in Fig [5|a. With finite- 
thickness bubble walls, however, it is necessary to let the z coordinates of points A and 
B approach ±00, and similarly to take the p coordinate of points D and E to infinity. 
Proceeding as in Ref. [1_1|, we realise that the line integral of eA^ is equal to the line 
integral of D^O, since the line integral of 8^9 along a closed path vanishes. Moreover, 
D a 6 must vanish on the line segments BD and EA, because they are far removed from 
the causal future of bubble intersections I + (C)t- Note, however, that D z 6 is non-zero on 
the segment ED. This is because the initial condition establishes a phase difference of 
26*o between D and E, which contributes to d z 6 a . In the axial gauge A z = 0, and we thus 
obtain 



DE/ 



- [2e - A6(t) 

e v 



(78) 



/AB 

where A# is given by Eq. ( [76|) or Eq. (|77|) . As a function of time, the Abelian flux 
rises from zero and performs damped oscillations about the late-time limit 29$/ e. This 



asymptotic value agrees with that found for the Kibble- Vilenkin solution in Sec. 3.1 
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7 Conclusions 



In this paper, we have attempted to provide the most complete analysis to date of the 
phase equilibration and field-strength evolution that occur in a U(l) gauge theory when 
two bubbles of true vacuum collide. Our analysis has been motivated by the work of 
Kibble and Vilenkin [11J , and we have been able to go beyond some of their assumptions 
concerning the initial configuration of the phase of the scalar field. In particular, we have 
found that a linear superposition ansatz, coupled with smooth bubble profiles arising 
from bounce solutions to the Euclidean equations of motion, leads to analytical solutions 
that are in excellent agreement with the full numerical evolution of the fields, at least for 
small initial phase difference between the two bubbles. 

In Sec. 4, we presented the new, smooth analytical solutions for the Abelian field 
strength and the gauge-invariant phase. They have a number of advantages over solutions 
obtained in Ref. [11| using the step-function approximation. For example, from the latter 
one obtains a peak value of the field strength that differs by a factor of four from the 
correct value obtained in numerical simulations (a factor of sixteen difference in energy 
density). Moreover, the solutions ( P^) - ([3H|) derived from the step-function approximation 
have discontinuities on the future null surface of bubble intersections and so are unable 
to account for the behaviour of the fields near that surface. The solutions are only 
valid over limited ranges of the coordinate r = \ft 2 — x 2 — y 2 , which means that they 
provide no information about the field dynamics shortly before the bubble collisions or, 
at later times, outside the bubble walls. The smooth ansatz that we adopt, on the 
other hand, allows us to account for the fact that realistic bubble walls of finite width 
have exponential tails that begin to overlap long before the instant of collision. The 
overlapping tails initiate the evolution of the phase of the scalar field and of the vector 
potential, so that it is already in an advanced stage when the collision reaches full impact. 
The advantage of the smooth approximation becomes manifest in the comparative plots 
of Figs. @-g. 

An interesting feature emerges in Fig. 5. There appears to be evidence for the ex- 
istence of alternating positive peaks and negative troughs in the field strength, corre- 
sponding to expanding circular flux-tube loops with alternating direction of the field 
strength. This raises the possibility that we are seeing here an analogy to the suite of 
alternating vortices and antivortices found by Digal and Srivastava [^9j in the collision 
region of two colliding bubbles in a global U(l) model. Our peaks are not large enough 
to form defects, because we have set the initial conditions so that the geodesic rule is 
obeyed, but it seems that our positive and negative peaks may be an "embryo" for such 
vortex-antivortex configurations. This is supported by the observation that the modulus 
of the Higgs field appears to have dips at the locations of the big peak and the trailing 
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anti-peak. In a hard collision, the field strength would be larger and the Higgs modulus 
would go down to zero in these two places, corresponding to the formation of topological 
defects. 

We should discuss some of the assumptions we have had to make in this article. 
We have restricted ourselves to 0(l,2)-symmetric solutions, but, as we have argued in 
Sec. |27^, this symmetry provides a good description in general for all pairs of nucleated 
bubbles except possibly in the very first stages of their expansion. Throughout our 
analysis we have kept the radial component of the complex scalar field fixed in the 
minimum of the potential within the union of the interior of the two bubbles. The validity 
of this approximation has been checked against the numerical simulations (see Fig. |j) 
and has been found satisfactory for the range of initial conditions we have adopted. The 
constraint placed upon us is that we are not really in a position to discuss the collision 
process involving large initial phase differences, for that is a situation where we would 
expect large gradient terms to arise and the scalar field to fluctuate wildly away from the 
vacuum. At the moment, we do not see any obvious way to approach this problem. 

The work we have presented here has a number of natural extensions. For exam- 
ple, one can use our approach to obtain smooth and realistic solutions to effective field 
equations that incorporate the effects of plasma conductivity, and thereby improve on 
previous results of this kind obtained with the step-function approximation ||T4|, O]. Al- 



ternatively, or perhaps in addition, one can try to include the effects of frictional damping 



of the bubble walls |21], [22|, 23"| . This unfortunately makes the problem non-relativistic 
and reduces the symmetry from 0(1,2) to 0(2). Therefore, it is likely that a new type 
of solution, or even approximation, will be needed. 

By going beyond the electroweak Standard Model, we could begin discussing the gen- 
eration of primordial magnetic fields in a first-order transition. We have good estimates 
for the magnitude and width of the magnetic fields formed from our smooth ansatz. This 
could be useful in determining the possible success of generating galactic-strength fields 
from such seeds in the early Universe. 
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Appendix A: Lattice Methods 



Here we discuss the method used for the numerical evolution of the equations of motion. 
The method employs the techniques of lattice gauge theory adapted to a Minkowski 



metric rather than a Euclidean metric. The formalism has been described by Moore |30 
for an SU(2) gauge group and was modified to suit our case of interest, a U(l) gauge 
theory. 

We approximate space-time by a Cartesian grid whose nodes are taken to have spatial 
separation a. Time is discretised with a fundamental time step ha, h < 1. Let us also 
define to be the lattice spacing in the \x direction, so that Iq = ha and U = a, i = 1, 2, 3. 
The quantity I is not a Lorentz vector, and so raising or lowering of indices is purely 
notational: Z M = Z M . A vector of length in the \x direction is denoted by fx. 

We proceed by introducing a unit-modulus complex vector field V : (//, x) i— > U(l) 
such that Vfj,{x) lives on the link between the neighbouring lattice sites x and x + fx. A 
scalar field <fi : x i— > C, with 4>(x) living on the node at x, completes the matter content of 
the model. The vector field living on the links between the sites alters the way the scalar 
field is transported around the lattice in much the same way as a connection in general 
relativity affects the transport of tensor quantities. The lattice covariant derivative of 
4>(x), Dfj,(f>, is defined by 



(d^){x) = i iv»(x)(j>(x + A) - . (A.i; 



Here there is no summation over repeated indices; in the following any summations shall 
be made explicit. The covariant derivative (|A. 1|) transforms in the same way as the scalar 
field under a local U(l) transformation provided that we have 

(j){x) — ► n{x)(j){x) , n{x) G U(l) ; 

V^x) — > Q(x)V^x)Q\x + fi) . (A.2) 

We then have the following candidate for a gauge-invariant kinetic term of the scalar 
field, 

^(d^ (x) (£>»<!>) (x) . (A.3) 

A gauge-invariant kinetic term for the complex vector field V^(x) can be found by realising 
that the quantity 

V^(x) = V v {x)V,{x + v)V}{x + fi)Vl{x) (A.4) 

is gauge- invariant. The relative normalisation of the scalar and vector kinetic terms 
is determined from the continuum limit a —>■ by making the identification V^{x) = 
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exp(— il^eAfj). In the limit of small a one obtains^ 



(D^)(x) = <9 M 0O) - ieA^x^x) + 0(a 2 



ReV^(x) 



1 - -IM ^ F ^ F ^ + ©(a ) > 
d/j,A u — d v A^ . 



(A.5) 

(A.6) 
(A.7) 



In terms of these quantities we may then write down a lattice action with a discrete 
version the gauge symmetry that reduces to the continuum U(l) gauge theory as a —>■ 0. 



S 



E ^ £ 



sites 
and links 



i (^(x + /^(x) - <p\x)) (V^xMx + A) - 0(x)) 



2 V- 1 



V„(x)V„(x + i>)Vj(x + jl)V}(x)] -V(\<f>(x)\) . (A.8) 



Here V is the Higgs potential for the scalar field 0. The equations of motion are derived 
by requiring that the lattice analogue of functional differentiation of S vanish with respect 
to all the fields. These functional derivatives can be calculated using the relations 



S(j)(x) 

SV^x) 
5V u (y) 

6V„{y) 



5 



x,y 



ha 4 



»ha 4 



-Vt(x)V}(x)5; 



ha A 



(A.9) 
(A.10) 
(A.ll) 



Upon doing so one discovers that the evolution of <p(x) is not fully determined. This is 
a consequence of the gauge-invariance of the system and can be remedied by making a 
gauge choice analogous to the temporal gauge, V = 1. The equations of motion are then 
found to be 



ir(x, t + At/2) = it{x, t - At/2) 
+ ha | E [ V i& t)<j>(x + i,t)- 2(f)(x, t) + V^(x - t, t)(f)(x - i, t) 



dV 



(A.12) 



6 In this Appendix we use a sign convention for the covariant derivative which is the opposite from 
that used in the rest of the paper. 
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lm[E h (x,t + At/2)) = hn[E k (x,t- At/2)) (A. 13) 

+ ha [— Im U>\x + k,t)V k \x,t)<f)(x,t) 

-^3 e (Vk(x, mix + k, mi{x + 1, m?(x, t) 

+ - i, t)V k (x, t)Vf(x + k-i, t)Vt(x - i, tj) 

In writing this we have found it convenient to define the quantities ir(x) and Ei(x) through 
the relations 

(j)(x,t + At) = <j)(x,t) + han(x,t + At/2) , (A.14) 
Vi(x,t + At) = ea 2 hE i (x,t + At/2)Vi(x,t) . (A.15) 



Eqns. ( |A.12| )-( [A.15| ) define a leap-frog algorithm where the canonical momenta, 7r(x) 



and Ei(x), are defined on time slices half-way between the time slices where the fields 
live. Note that ( |A.13| ) only shows how to evolve the imaginary part of Ei. To get 
the full evolution we use the fact that V{ has unit modulus; Eq. ( |A.15| ) therefore gives 



\Ei\ = l/(ea 2 h), which fixes the real part of Ei up to a sign. To determine the sign we 
write Ek = Ef c + iE l k m and substitute V k = exp(—ieaA k ) into (|A.15|) , finding to lowest 
order in h: 

exp (-ieaA k - iea 2 hA k ^ ~ ea 2 h (Ef c + iE l k m ^) exp (-ieaA k ) , (A. 16) 

i.e., 

1 - iea 2 hA k ~ ea 2 h (Ef c + iE l k m ) . (A.17) 



We find therefore that in the small-/;, limit the sensible sign for the real part of eE k is 
positive. The sign of the charge e is defined by matching Eq. (|A.5| ) to the usual covariant 



derivative in the continuum limit. This completes the evolution equations. The equation 
found by varying the action with respect to Vq gives the following lattice Gauss constraint, 

^E Im i E i( x > 1 + At / 2 ) - Ei{x -i,t + At/2)} = 



2 elm 



7r\x,t + At/2)(f>(x,tj\ , (A.18) 



which is naturally satisfied because the action is invariant under the discrete version of 
the gauge symmetry, Eq. ( |A.2|) . 

When working on the lattice one must find appropriate lattice quantities which reduce 
to the objects under study in the continuum limit. In this work we have been particularly 
interested in the variation of the phase of the scalar field. The gauge-invariant gradient 
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of the phase is essentially equal to the Nother current. This gives us the motivation for 
defining the lattice covariant derivative of the phase 

be = ^{{b<p)^-^b<p) , (A.i9) 

which is clearly gauge-invariant. For the lattice version of the Abelian field strength we 
use the relation 

ImVr = el^F^ + 0(a A ), (A.20) 
which motivates the lattice definition of the field strength as 

Bi = -^—e ijk lmV jk . (A.21) 



Appendix B: Poly logarithms 

The poly-logarithm function Li n (z) is defined by 

(-l) n 



LUz) 



1 , ln n_1 x 
dx- 



oo z k 



(n — 1)! jo x — 1/z 



E 

k=l 



(b.i; 



For real z and integer n the power series is convergent for — 1 < z < 1, n > and for 
z — 1, n > 1. The functions Li n (z) are real- valued for z < 1. In the case n = 2 it is 
convenient to use the functional relations 



U 2 (z) 



71 



Li 2 (l - z ) + — - lnzln(l - z) 
o 

1 - z 



Li 2 
- Li 2 



7T 1 o 

, h - hr 

i-^y 6 2 |- 



- -ln 2 H 



±<z<l 

2 < -1 

-1< ^ < 



(B.2) 



to express Li 2 (z) in terms of functions Li 2 (w) where < w < 1/2. For such w the power 
series in Eq. ( B.I ) converges rapidly. For general complex z we also have the following 
relations: 

(B.3) 

IiaC*) + Li 2 (-) = - i ln 2 (-z) , (B.4) 



Li 2 (z) + Li 2 (l — z) = lnzlnfl — z) 

6 
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From Eqs. ( [B.4j ) and ( |B . 1| ) it follows that the dilogarithm function Li2 has the following 
asymptotic behaviour: 

1 7T 2 

Li 2 (x) -> -- ln 2 (-x) - — , ar -> -oo . (B.5) 



More relations between Li 2 functions are given in Ref. |28| . 

References 

[1] A. Vilenkin and E.P.S. Shellard, Cosmic Strings and Other Topological Defects 
(Cambridge University Press, 1994). 

[2] T. Vachaspati and G. E. Volovik, Int. J. Mod. Phys. B 10, 471 (1996). 

[3] A. Guth, Phys. Rev. D 23, 347 (1981). 

[4] Daile La and P. Steinhardt, Phys. Rev. Lett. 62, 376 (1989); Erratum ibid. 62 1066 
(1989). 

[5] D. Lyth and E. Stewart, Phys. Rev. Lett. 75, 201 (1995). 

[6] Bao-Hua Liu, L. McLerran, and N. Turok, Phys. Rev. D 46, 2668 (1992). 

[7] M. Laine and K. Rummukainen, Phys. Rev. Lett. 80, 5259 (1998); J.M. Cline and 
G.D. Moore, Phys. Rev. Lett. 81, 3315 (1998). 

[8] S. Hawking, I. Moss, and J. Stewart, Phys. Rev. D 26, 2681 (1982). 

[9] S. Rudaz and A. Srivastava, Mod. Phys. Lett. A 8, 1443 (1993). 

[10] M. Hindmarsh, A.-C. Davis, and R. Brandenberger, Phys. Rev. D 49, 1944 (1994). 

[11] T. Kibble and A. Vilenkin, Phys. Rev. D 52, 679 (1995). 

[12] E. Copeland and P. Saffin, Phys. Rev. D 54, 6088 (1996). 

[13] P. Saffin and E. Copeland, Phys. Rev. D 56, 1215 (1997). 

[14] T. Kibble, J. Phys. A 9, 1387 (1976). 

[15] I. Chuang, R. Durrer, N. Turok, and B. Yurke, Science 251, 1336 (1991). 

[16] J. Ye and R. Brandenberger, Mod. Phys. Lett. A 5, 157 (1990); A. Srivastava, Phys. 
Rev. D 45, R3304 (1992); ibid. D 46, 1353 (1992). 

40 



[17] K. J. Ahonen and K. Enqvist, Phys. Rev. D 57, 664 (1998). 
[18] O. T6rnkvist, Phys. Rev. D 58, 043501 (1998). 



[19] O. Tornkvist, [hep-ph / 9902432| , in Proceedings of the International Conference on 



Strong and Electroweak Matter (SEWM-98), Copenhagen, Denmark, 2-5 Dec 1998 
(World Scientific). 



[20] A.-C. Davis, M. Lilley, and O. Tornkvist, Phys. Rev. D 60, 021301 (1999), |astro- 



ph/9904022 . 



[21] A. Ferrera and A. Melfo, Phys. Rev. D 53, 6852 (1996). 

[22] M. Joyce, T. Prokopec, and N. Turok, Phys. Rev. D 53, 2930 (1996). 

[23] A.-C. Davis and M. Lilley, Damped Bubble Collisions in First-Order Phase Transi- 
tions, Cambridge preprint DAMTP-1999-87. 

[24] S. Coleman, Phys. Rev. D 15, 2929 (1977); Erratum ibid. D 16, 1248 (1977). 

[25] J. Garriga, Phys. Rev. D 49, 5497 (1994). 

[26] B. Gidas, W. Ni, and L. Niremberg, Comm. Math. Phys. 68, 209 (1979). 

[27] B. de Wit and J. Smith, Field Theory in Particle Physics, Vol. 1 (Elsevier, Amster- 
dam, 1986), Appendix D. 

[28] L. Lewin, Polylogarithms and Associated functions (North Holland, 1981). 

[29] S. Digal and A.M. Srivastava, Phys. Rev. Lett. 76, 583 (1996). 

[30] G. D. Moore, Nucl. Phys. B 480, 657 (1996); ibid. B 480, 689 (1996). 

[31] H. Kurki-Suonio and M. Laine, Phys. Rev. Lett. 77, 3951 (1996). 



41 



